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Abstract 

Background: Yersinia enterocolitica outer membrane protein A (OmpA) is one of tlie major outer membrane 
proteins witli liigli immunogenicity. We performed tlie polymorpliism analysis for tlie outer membrane protein A 
and putative outer membrane protein A [p-ompA] family protein gene of 318 Y. enterocolitica strains. 

Results: The data showed all the pathogenic strains and biotype lA strains harboring ystB gene carried both ompA 
and p-ompA genes; parts of the biotype lA strains not harboring ystB gene carried either ompA or p-ompA gene. In 
non-pathogenic strains (biotype lA), distribution of the two genes and ystB were highly correlated, showing genetic 
polymorphism. The pathogenic and non-pathogenic, highly and weakly pathogenic strains were divided into different 
groups based on sequence analysis of two genes. Although the variations of the sequences, the translated proteins 
and predicted secondary or tertiary structures of OmpA and P-OmpA were similar. 

Conclusions: OmpA and p-ompA gene were highly conserved for pathogenic Y. enterocolitica. The distributions of two 
genes were correlated with ystB for biotype 1 A strains. The polymorphism analysis results of the two genes probably due 
to different bio-serotypes of the strains, and reflected the dissemination of different bio-serotype clones of Y. enterocolitica. 
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Background 

Y. enterocolitica is spread primarily through contami- 
nated food or water [1] and causes a wide range of intes- 
tinal diseases, including enteritis, mesenteric lymphadenitis, 
and sepsis in some severe cases; and also cause some 
complications such as erythema nodosum and reactive 
arthritis [2]. Currently, Y, enterocolitica is divided into 
six biotypes (lA, IB, and 2-5) and more than 50 sero- 
types [3]. The biotypes of Y, enterocolitica are divided 
into three groups according to the bacterial pathogenic 
properties: non-pathogenic biotype lA, weakly patho- 
genic biotypes 2-5, and highly pathogenic biotype IB 
[4]. At present, the virulent factors of Y, enterocolitica 
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are mainly referred to type III secretion system (TTSS) 
encoded by 70-kb plasmid pYV [3,5], Yersinia adhesin 
A (YadA), the virulence genes involved inv (invasion gene), 
ail (attachment and invasion locus gene), yst {Yersinia 
stable toxin gene), myfA (mucoid Yersinia factor gene) and 
the pathogenic island (HPI), etc. [6]. Biotype lA strains are 
traditionally considered non-pathogenic, however in recent 
studies have confirmed that a portion of them can cause 
clinical symptoms similar to pathogenic strains [7]. 

Bacterial outer membrane proteins primarily contain 
outer membrane protein A, porin C and F, being the 
major immunogenic proteins, and widely present in 
Gram-negative Enterobacteriaceae [8]. Recently, multi- 
function of outer membrane protein A of intestinal 
bacteria has been demonstrated [9-12], but seldom re- 
ferred to Y, enterocolitica. In our previous study [13], 
OmpA was the major immunogenic protein of both 
highly and weakly pathogenic Y, enterocolitica incu- 
bated at different temperatures. To further identify the 
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characteristics of ompA for Y. enter ocolitica, we se- 
quenced and analyzed the polymorphism of ompA (and 
P'OmpA) genes of Y, enter ocolitica. 

Results 

Distribution of ompA and p-ompA genes 

The ompA and p-ompA genes were detected in 318 
Y. enterocolitica strains. The data showed 170 of all the 
pathogenic strains carried both ompA and p-ompA; 91 
biotype lA strains carried ompA, and 106 carried p-ompA 
(Table 1). The ompA and p-ompA genes were both existed 
for biotype 1 A strains carried ystB gene. However, the iso- 
lates only had one of the two genes for biotype lA strains 
without ystB gene. The distribution of the two genes was 
associated with ystB, especially for ompA (Tables 2 and 3). 
The correlation analysis (P < 0.05) showed correlation 
coefficient (r) was 0.67 between ompA and ystB and 
0.58 between p-ompA and ystB, 

Polymorphism analysis of the two genes 

ompA: The Open Reading Frame (ORF) of ompA was 
1,074 bp encoding 357 amino acids (Genbank: YP_00 
1005874.1). The ompA gene of 261 strains formed 23 se- 
quence types. The pathogenic group contained five 
types. Pattern A-E (Figure lA). 155 pathogenic strains 
(76 bio-serotype 3/0: 3, two 4/0: 3, 68 2/0: 9, one 4/0: 
9 and eight 3/0: 9) were clustered into Pattern A, and 
accounted for 91.2% of all of the pathogenic strains. 
Therefore, pattern A was an absolute primary type in all 
pathogenic isolates. Other nine pathogenic strains were 
clustered into pattern B and C, seven (two bio-serotype 
2/0: 3, four 4/0: 3 and one 3/0: 9) were clustered into 
Pattern C, two pathogenic bio-serotype 2/0: 5, 27 clus- 
tered into pattern B. Few nucleotide differences were 
found between pattern A and B (Figure 2); while the nu- 
cleotide insertion was found compared with pattern A and 
C (Figure 3). Six highly pathogenic bio-serotype IB/O: 8 
strains were clustered into Patterns D and E (Figure lA), 
and the nucleotide differences were shown with red bases 
compared with pattern A and B (Figure 2). 

The non-pathogenic group contained 91 biotype lA 
strains, formed 18 sequence types. Pattern F-W; all the 83 

Table 1 Distribution of ompA and p-ompA gene in 
all strains 



ompA+ ompA- Total 





p-ompA+ 


p-ompA- 


p-ompA+ 


p-ompA- 




Pathogenic strains 


170 


0 


0 


0 


170 


Biotype lA 
strains ystB+ 


84 


0 


0 


0 


84 


Biotype lA 
strains ystB- 


0 


7 


22 


35 


64 


Total 


254 


7 


22 


35 


318 



+: positive; -: negative. 



Table 2 Distribution of ompA and ystB gene in biotype 
1A strains 



ompA 


+ 


ystB 




Total 


+ 


84 




7 


91 




0 




57 


57 


Total 


84 




64 


148 



+: positive; -: negative. 



strains carried the ystB gene were clustered into pattern 
F-T; Eight strains (except one) without ystB gene formed 
pattern U-W (Figure lA). 

There were nucleotide fragment insertions and dele- 
tions for ompA gene sequences. The nucleotide fragment 
insertion was found for pattern C. The insertion site lo- 
cated at 659-664 nt of the ORF with the "-CTCCAG-" 
compared with pattern A (Figure 31), this made an ala- 
nine (A) and a proline (P) inserted at position 220 and 
221 in the amino acid sequence of the OmpA. For bio- 
type lA strains, sequence type diversities were found 
compared with pathogenic strains, and the nucleotide 
fragment deletions were observed for pattern U to W. 
15 nucleotides deletions located at 140-154 nt with 
"-TCTATGGCGACCAAG-", and 6 nucleotides deletions 
located at 409-414 nt with "-GACGCA-" of ompA were 
found (Figure 3D). Eventually, this led to a change of the 
amino acid of the OmpA at sites 39-52, 94-98,133-146, 
included amino acid deletions or translation into other 
amino acids. 

Although some predicted amino acid differences were 
found for ompA of Y enterocolitica, the primary proteins 
and its functions were the same, even if the insertions of 
the pathogenic or deletions of the non-pathogenic strains 
occurred. For all the Y. enterocolitica possessed ompA in 
this study, the translated proteins were identical. The pre- 
dicted secondary structure of OmpA for all the strains 
were almost the same, and the predicted tertiary structure 
of OmpA for all the isolates were similar as well. It was 
highly conserved for the structure and function of OmpA 
for Y, enterocolitica, no matter the pathogenic ability and 
other characteristics. 

P'OmpA: The Open Reading Frame (ORF) of p-ompA 
was 1,377 bp encoding 458 amino acids (GenBank: 



Table 3 Distribution of p-ompA and ystB gene in biotype 
1A strains 



p-ompA 


+ 


ystB 




Total 


+ 


84 




22 


106 




0 




42 


42 


Total 


84 




64 


148 



+: positive; -: negative. 
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PatternCode / No, 
^ Pattern W / 3 



Pattern V / 2 


1 Pattern U /3| 


Pattern 


T 


/13 ■ 


Pattern 


S 


/I 


Pattern 


R 


/3 


Pattern 


Q 


12 


Pattern 


P 


/I 


Pattern 


0 


Ml 


Pattern 


N 


/11 


Pattern 


M 


/1 


Pattern 


L 


/1 


Pattern 


K 


/I . 


Pattern 


J 


/3 " 


Pattern 


1 


/1 _ 


Pattern 


H 


/4 ■ 


Pattern 


G 


/I 


Pattern 


F 


/1 



Non-pathogenic Group 
(Biotype ^^lystB-) 

K 

I Containing 1 Biotype WystB+ | 



OOj- • Pattern E / 4 



100| ^# 



Pattern D / 2 _ 

• Pattern C / 7 " 

Pattern B 12 

Pattern A /155 



Non-pathogenic Group 
(Biotype ^NystB+) 



Non-pathogenic Group 
(Biotype ^A^ystB+) 

Non-pathogenic Group 
(Biotype ^f\lystB+) 

Pathogenic Group 
(Biotype 1B) 

Pathogenic Group 
(Biotype 2-4) 



99j£ 



100 



lOOj- 

rt 

k 



4 



PatternCode / No. 

Pattern 38/1 ' 
Pattern 37 / 1 
Pattern 36 / 6 
Pattern 35 / 2 
Pattern 34 / 1 
Pattern 33 / 1 
Pattern 32 / 1 
Pattern 31 / 2 
Pattern 30 / 1 
Pattern 29 / 1 
Pattern 28/ 2 . 
Pattern 27 / 5 ' 
Pattern 26 / 7 
u Pattern 25 / 1 
t Pattern 24 / 1 
' Pattern 23 / 1 
Pattern 22 / 2 
© Pattern 21 / 4 
0 Pattern 20 / 3 
Pattern 19/ 29 
Pattern 18/1 
Pattern 17/3 
Pattern 16/1 
Pattern 15 /I 
Pattern 14 / 1 
Pattern 13/1 
Pattern 12/9 
Pattern 11/3 
Pattern 10/9 
Pattern 9 / 1 
Pattern 8 / 1 . 

• Pattern 7/81' 
Pattern 6 / 1 
Pattern 5 / 4 

• Pattern 4 / 2 



Pattern 3 / 79 



_y«|-# Pattern 
Pattern 



2 /5 
1 /1 

B 



Non-pathogenic Group 
(Biotype ^NystB-) 



Non-pathogenic Group 
(Biotype ^AIystB+) 



Pathogenic Group 
(Biotype 1 B,2-4) 



J (Biotyi 
I Containing 3 Bio-serotype 1A/0:3T 



Figure 1 Cluster tree of ompA and p-ompA gene sequences. A: Cluster tree of ompA gene sequences from 261 strains; B: Cluster tree of p-ompA 
gene sequences from 275 strains; red: pathogenic strains; green: non-pathogenic strains. 
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Figure 2 Sequence polymorphisms of ompA gene for pathogenic strains. The number above bases represented the position of bases in the ORF; 


figure in brackets represented strain number; red represented mutant bases; yellow area represented sense mutations, and others were nonsense mutations. 
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Figure 3 ompA nucleotide insertions and deletions. N: the ORF of pattern A for ompA; D: the deletion ORF of pattern U-W for ompA; I: the insertion 
ORF of pattern C for ompA; blue: nucleotide differences concentrated area; yellow: nucleotide deletion position; green: nucleotide insertion position. 



YP_001006877.1). P-ompA genes were existed among 276 
Y. enter ocolitica, 170 pathogenic and 106 biotype lA 
strains (84 carried ystB gene and 22 without). The p-ompA 
genes were clustered into 38 sequence types, divided into 
pathogenic and non-pathogenic group (Figure IB). Pattern 
1-7 belonged to pathogenic group, included all the 170 
pathogenic strains and three bio-serotype lA/O: 3 isolates. 
Pattern 3 and pattern 7 were the primary types for the 



pathogenic strains, 67 bio-serotype 2/0: 9, eight 3/0: 9, 
one 4/0: 9 and three bio-serotype lA/O: 3 strains formed 
pattern 3; 75 bio-serotype 3/0: 3, one 2/0: 3, four 4/0: 3 
and one 3/0: 9 formed pattern 7. Pattern 5 contained two 
bio-serotype 4/0: 3, one 3/0: 3 and one 2/0: 9 strains; 
pattern 6 contained only one 2/0: 3 isolate; two bio- 
serotype 2/0: 5, 27 strains formed pattern 4, and six 
highly pathogenic IB/O: 8 strains formed pattern 1 and 2. 
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Figure 4 Sequence polymorphisms of p-ompA gene for pathogenic strains. The number above the bases represented position of the bases 
in the ORF; figure in brackets represents the strain number; red represented mutant bases; yellow area represented sense mutations, and others 
were nonsense mutations. 
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Some nucleotide point mutants were found among pattern 
1 to 7, as shown in Figure 4. 

Pattern 8-38 referred to non-pathogenic group, in- 
cluded 106 biotype lA strains. All of the 84 isolates carried 
ystB were clustered into pattern 8 to 27; 22 strains without 
ystB were clustered into pattern 28 to 38 (Figure IB). 

Although the diversity of cluster results was found for 
P'OmpA, the primary proteins structure and function of 
different patterns were identical predicted by software. 
The predicted secondary and tertiary structure or func- 
tions of proteins for all patterns were also similar, shown 
the conservative characteristic. 

Discussions 

The surface structures of bacterial pathogens (including 
pilus, flagella, outer membrane proteins, and various se- 
cretion systems) are likely to interact with host tissue to 
regulate adhesion and invasion [14]. The outer membrane 
protein A belong to highly conserved protein in intestinal 
bacteria, and play a key role in bacterial integrity and 
virulence [15]. Currently, more evidence shows the 
pathogenicity of outer membrane protein A in a variety 
of pathogens [11,16-19]. 

In our study, all the 170 pathogenic Y, enterocolitica 
strains carried the ompA and p-ompA genes. Few nucleo- 
tides changes were found for both ompA and p-ompA of 
the pathogenic strains. Therefore, the distributions of two 
genes were highly conserved for pathogenic Y. enterocolitica. 
The translated proteins or predicted structures of dif- 
ferent patterns of the two genes were the same, certifi- 
cated the conservative property of ompA and p-ompA 
for Y. enterocolitica. Several researches had been widely 
shown the distribution of outer membrane protein A 
among entero-pathogenic bacteria, and its important 
role in bacterial infection and immunogenicity [16-19]. 
However, seldom study referred to ompA or p-ompA of 
Y, enterocolitica, so it was the first time to perform this 
research. Our results showed the pathogenic strains 
and non-pathogenic strains were divided into different 
groups, and highly or weakly pathogenic strains were 
also distinguished based on sequence results of two 
genes, which reflected the different bio-serotype distribu- 
tions of Y, enterocolitica. In China, serotype O: 3 and O: 9 
strains were predominate pathogenic Y, enterocolitica y and 
most of these bacteria isolated from patients, swine and 
dogs. While, non-pathogenic strains referred to biotype 
lA, and widely distributed among different hosts. Further- 
more, no highly pathogenic IB/O: 8 strain was isolated in 
China up to present, all the highly pathogenic IB/O: 8 
were foreign strains in our study. Additionally, the weakly 
pathogenic strains (biotype 2-4, serotype O: 3 or O: 9 
strains) included wild strains from China and reference 
strains not from China showed no differences even if they 
were isolated from different origins in our study. Therefore, 



the cluster results for ompA or p-ompA explained the dif- 
ferent bio-serotype distributions of Y, enterocolitica. Highly 
pathogenic biotype IB strains have been shown to differ 
genetically from weakly pathogenic biotype 2-4 strains, 
and they belonged to different subtypes. The poly- 
morphism analysis results of the two genes probably 
reflected the dissemination of different bio-serotype 
clones of Y, enterocolitica for a period of time. 

Biotype lA strains lack pYV plasmid and chromosomal 
virulence genes, and generally regarded as avirulent [1]. 
However, few studies have conflrmed biotype lA strains 
were related to outbreaks of nosocomial infections and 
foodborne diarrhea [20,21]; and some early studies found 
that biotype lA strains could cause abortion in goats 
and cattle [22-25]. Grant et al. [26,27] showed biotype 
lA strains invaded epithelial cells and resisted the kill- 
ing effect of macrophage. Also biotype lA strains were 
associated with the potential pathogenicity in humans 
[28]. Enterotoxin is an important pathogenic factor in 
most enteric pathogens, and the ystB gene coded for a 
class of thermo-resistant enterotoxin in biotype lA Y. 
enterocolitica [29,30]. Virulence related gene ystB was a 
distinguishing marker of biotype lA strains, presented 
in close to 100% of clinical isolates [31,32]. However, 
whether ystB gene as a virulent factor for biotype lA of 
Y, enterocolitica has not been confirmed. Nakano et al. 
[33] found Salmonella enterotoxin (stn) regulated the 
OmpA membrane localization and functions, indicated 
the close relationship between them. In our research, 
ompA and p-ompA were correlated with ystB in biotype 
lA strains, and formed the independent cluster patterns, 
the strains with ystB or without ystB for biotype lA were 
also separated, which indicated the phenomenon that 
OmpA was linked with enterotoxins for Y, enterocolitica 
biotype lA strains. 

Conclusions 

We showed the ompA and p-ompA genes of Y, enterocolitica 
were highly conserved in pathogenic strains; specially, 
the two genes showed a high correlation with ystB in 
biotype lA strains. The pathogenic and non-pathogenic 
strains, highly and weakly pathogenic strains were divided 
into different groups based on sequence polymorphism 
analysis of the two genes, which reflected the different 
bio-serotype distributions of Y, enterocolitica. 

Methods 

Bacterial strains and identification of biotype and serotype 

The bacterial strains used in this study were screened 
from the Chinese Yersinia enterocolitica library which 
contains nearly 4,000 strains gathered by our laboratory 
and derived from diarrhea patients, food, animals, and the 
environment. Strains were selected to cover different iso- 
lation dates, different hosts, and separated locations. We 



Li et a I. BMC Genomics 2014, 15:201 
http://www.bionnedcentral.conn/1471 -21 64/1 5/201 



Page 6 of 8 



Table 4 The information of Y, enterocolitica used in this study 



Source 



Pathogenic strains (bio-serotype) 



Non-pathogenic strains (bio-serotype) 
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a: K enterocolitica W22703, contig 7180000001374, GenBank: FR71 8562.1; b: Y. enterocolitica subsp. palearctica 105.5R (r), complete genome, GenBank: 
CP002246.1; c: Y. enterocolitica subsp. palearctica Y11, GenBank: FR729477.2; d: Y. enterocolitica subsp. enterocolitica 8081 complete genome, GenBank: 
AM286415.1; UN: undetermined serotype. 



chose 150 pathogenic and 148 biotype lA Y. enterocolitica 
strains isolated from China; 16 pathogenic reference 
strains from Europe, United States, and Japan; and four 
pathogenic complete-genome-sequenced strains (Table 4). 
The serotypes of these strains were determined as previ- 
ously described [1,34-36], and the biotypes of strains were 
identified using the scheme reviewed by Bottone [37] . The 
pathogenic strains were positive for all genes {ail^, ystA^, 
virF^y and yadA^); however, some pathogenic strains lost 
the plasmid virulence genes for virF and yadA, but still 
had ail and ystA genes located on the chromosome, the 
non-pathogenic strain was negative for all these genes. 

The sample collection and detection protocols were ap- 
proved by the Ethics Review Committee from the National 
Institute for Communicable Disease Control and Preven- 
tion, Chinese Center for Disease Control and Prevention. 



Primer design 

Two genes for OmpA of Y. enterocolitica reference strain 
8081 (NC_008800.1) were shown when we searched the 
NCBI web, one was ompA, another was p-ompA, There- 
fore, we designed the primers of the two genes by using 
CloneManager software 4.0, and the primers were showed 
in Table 5. Primers were synthesized by Shanghai Sangon 
Biological Engineering & Technology and Service Co., 
Ltd, China. 

PGR, DNA sequencing and sequence analysis 

Bacteria were cultured as previously described [35]. The 
bacterial DNA was extracted using a Blood & Tissue Kit 
(QIAGEN, USA). PCR was performed in a 20 [A volume 
containing 10 [A PCR premix (TaKaRa, Japan), 8 [x\ ultra- 
pure water, 0.5 [A of each forward and reverse primer 



Table 5 Primers and annealing temperatures for ompA and p-ompA 



Target gene 


Primer direction 


Primer Sequences (5' ^3') 


GenBank no. 


Location 


Amplicon lengt 


Annealing temp 


ompA 


Forward 


ACATCACAOTGTAACmCTCACC 


YP_001 005874.1 


1783285-1783261 


1451 bp 


58°C 




Reverse 


AGAAGTATCAGAATCAGATGTCGTC 




1781835-1781859 






p-ompA 


Forward 


GCGGCAAATOCGTACAGTG 


YP_001 006877.1 


2919405-2919386 


1560 bp 


60°C 




Reverse 


CAGCCCACCAGCAATATOG 




2917806-2917825 
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(25 (imol/l), and 10 ng DNA template. Thermal cycling 
was performed in a MJ PTC200 (Bio-Rad, USA) and the 
conditions were: denaturation at 94°C for 5 min, followed 
by 25 cycles of melting at 94°C for 25 s, annealing for 30 s 
at various temperatures depending on the primers used 
(Table 5), elongation at 72°C for 30 s, and a final extension 
at 72°C for 10 min. The specific PGR products were 
purified using a Gel Extraction Kit (QIAGEN, USA) 
and sequenced at TaKaRa Biotechnology (Dalian) Go., 
Ltd. Nucleotide sequence alignments and cluster tree 
construction were performed using MEGA (Version 4). 
The statistical tests were performed using statistical ana- 
lysis software SAS version 9.2 (Statistics Analysis System). 

The different sequences of two genes were translated to 
predict amino acid by MEGA 4.0 software, and the second 
structure of the proteins were predicted by PredictProtein 
(http://www.predictprotein.org); the tertiary structure of 
proteins were predicted and analyzed by SWISS -MODEL 
(http://swissmodel.expasy.org/workspace). 

Availability of supporting data 

All types of patterns for ompA and p-ompA of Yersinia 
enter ocolitica in our study were uploaded to LabArchives 
(http://www.labarchives.com/bmc) Electronic Laboratory 
Notebook. All sequences can be shared from the following 
links: 

ompA sequences: 
https://mynotebook.labarchives.com/share/hqjing/ 
MjAuOHwzMjk5OC8xNi0yLlRyZWVOb2RlLzMwN- 
jA0MzQxNDh8NTIuOA=DOI:10.6070/H4MP517C. 

P'OmpA sequences: 
https://mynotebook.labarchives.com/share/hqjing/ 
MjIuMXwzMjk5OC8xNy00LlRyZW- 
VOb2RlLzQwMDM0NDM2MTB8NTYuMQ= 
DOI:10.6070/H4GX48HN. 
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